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INJECTION OF AN ATTACHED INVISCID JET AT AN OBLIQUE 
ANGLE TO A MOVING STREAM 
by Marvin E. Goldstein and Willis Braun 
Lewis Research Center 

SUMMARY 

An analytical solution has been obtained to the problem of a two-dimensional invis- 
cid, incompressible jet injected into a moving stream from an orifice set at an oblique 
angle to the stream for the case where the jet does not separate from the downstream 
edge of the orifice. The solution is valid when the difference between the total pressure 
in the jet and the total pressure in the main stream is not too large. Typical flow pat- 
terns are presented to illustrate the effects of varying both the orifice angle and the total 
pressure within the jet. 


INTRODUCTION 

The flow field resulting from the oblique penetration of a jet into a flowing stream is 
of considerable interest in a number of fluid mechanical devices. Among these are 
ground- effects machines, jet flaps, wing fans, and fuel injection systems. 

The dynamics of jet injection into moving streams is by no means fully understood. 
However, a certain amount of insight into this phenomenon can be gained by considering 
the injection of two-dimensional inviscid jets into flowing streams since flows of this 
type are simple enough to be amenable to mathematical analysis. It is realized that 
viscous effects can be significant in real flows. In order to take viscous effects into ac- 
count, however, it is necessary to first perform an inviscid analysis and then modify the 
flow by superimposing viscous boundary layers. In any event, it is hoped that the invis- 
cid analysis will reveal some of the significant features of the flow and thereby lead to 
an increased understanding of the phenomena involved. The work done along these lines 
to date is summarized in reference 1. 

In reference 1 a technique was developed and applied to obtain an explicit analytic 
solution for the flow field resulting from a two-dimensional inviscid and incompressible 


jet issuing from an orifice into a moving stream. The orifice was at an oblique angle to 
the stream and it was assumed that the jet separates from the downstream edge of the 
orifice to form a stagnant wake. Since it is not possible to tell from the inviscid analy- 
sis whether the jet will separate, remain attached, or form a separation bubble, it is 
also of interest to consider the flow that would result if no separation occurred. Hence, 
in this report the techniques used in reference 1 are used to obtain a solution for the 
case where the jet does not separate from the downstream edge of the orifice. The flow 
configuration is shown in figure 1 (p. 5). Since no separation is allowed to occur, it is 
necessary to allow the velocity to become infinite at the downstream edge of the orifice. 
This condition is approximately realized in certain real flows. The usual method for 
handling this situation is to replace the boundary streamline by one which is close to it 
when interpreting the results (ref. 2). In the present case this corresponds to replacing 
the infinitely thin plate which forms the downstream edge of the orifice by one of finite 
thickness with a streamlined leading edge. 

As in reference 1, the flow will be assumed to be two-dimensional, inviscid, and in- 
compressible. In addition, it will be required that in a certain sense (to be specified 
more precisely) the difference between the total pressure in the jet and the total pres- 
sure in the main stream be small. The upstream boundary of the jet is the streamline 
emanating from the upstream edge of the orifice. 

The problem is solved by expanding the solutions in a small parameter related to 
the difference in total pressure between the jet and the mainstream. The zeroth- order 
solution corresponds to equal total pressures. 

Since the boundary shapes for the first-order problem are unknown, a technique 
similar to that employed in thin airfoil theory is used to transform the first-order bound- 
ary conditions to the zeroth-order boundary. The solution is then obtained by using the 
theory of sectionally analytic functions. 

It is proved in the appendix of reference 1 that the solution obtained herein also ap- 
plies to the case where the density of the fluid in the jet differs from that in the main- 
stream. 


SYMBOLS 

A horizontal distance between edges of orifice 
a A/l 

B vertical distance between edges of orifice 
b B/l 
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C pressure coefficient along slip line 

ir° 

D* flow regions in physical plane 
Qq regions in T-plane 

H asymptotic jet width 

h H /l 

I function defined by eq. (54) 

J function defined by eq. (55) 

l characteristic length (set equal to Hq) 

M function defined by eq. (63) 

O order symbol 

o order symbol 

P total pressure in main stream 

Pj total pressure in jet 

P^ total pressure in main stream 
P.V. Cauchy principal value 

p static pressure 

p Q static pressure at jet source (far inside the orifice) 

p TO static pressure far upstream from jet 

Q volume flow through jet 

S slip line in physical plane 

S distance along slip line 

S^q slip line in T-plane 

T intermediate variable, T = £ + if] 

TJ X- component of velocity 

u U/V TO 

V Y- component of velocity 

V* velocity along slip line inside of jet 
s 

free stream velocity 

v/v^ 


V 



W dimensionless complex potential, cp + 
X coordinate in physical plane 
x X/l 

Y coordinate in physical plane 

y y/i 


z dimensionless complex physical coordinate, x + iy 

s 

z dimensionless coordinate of points on slip line 
r function defined by eq. (49) 
y dummy variable to replace 77 



A location of downstream edge of orifice in T-plane 
5 defined in fig. 7 

Z dimensionless complex conjugate velocity, u - iv 
r\ coordinate in T-plane 
0 function defined by eq. (31) 

A defined by eq. (58) 

4 coordinate in T-plane 
p density 

r dummy variable in T-plane 
$ velocity potential 
cp &/IV" 

^ stream function 

* v/iv* 

Subscripts: 

0 zeroth-order quantity 

1 first-order quantity 
Superscripts: 

5 value of quantity on slip line 
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+ value of quantity inside jet and orifice 
value of quantity in main stream 
(overbar) complex conjugate 


ANALYSIS 

Formulation and Boundary Conditions 


It will be assumed that the flow is inviscid, incompressible, and irrotational. The 
jet configuration is illustrated in figure 1. The analysis is limited to the case in which 
the difference between the total pressure in the jet Pj and the total pressure in the main 
stream P w is not too large; or more specifically, to the case in which 

Id «» 


where 



a) 


p is the density of the fluid and is the velocity of the main stream at infinity. 



Figure 1. - Jet penetrating stream. 
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Let Z be a convenient reference length which will be specified in the course of the 
analysis. The X and Y components of the velocity, U and V, respectively, will be 
made dimensionless by V^; and the stream function Sk and the velocity potential $ 
will be made dimensionless by V m l. Thus, the dimensionless quantities u, v, xp, and 
cp are defined by 


u 



v = 



* = 




V oo l 


<P 


— <& 


VI 


The dimensionless complex conjugate velocity £ and the dimensionless complex poten- 
tial W are defined, as usual, by 


£ = u - iv 


and 


W = cp + ixp 


With all lengths made dimensionless by Z (i. e. , x = X/Z, y = Y/l , a = A/Z, b = B/Z, 
and h = H/Z) the flow configuration is shown in the physical plane (with the complex var- 
iable z defined by z = x + iy) in figure 2. 

The stream of fluid issuing from the orifice formed by the two parallel walls HD 
and EH meets the main stream at the point D and forms a common streamline which 
is denoted by S in figure 2. The jet does not separate from the wall EH at the point 
E but turns and flows back along EC. As a result of using this model it is necessary 
to allow the velocity to become infinite at the point E. Points on the common stream- 

s s s 

line will be denoted by z = x + iy . 

In order to satisfy the requirement that there be no discontinuities in static pressure 
anywhere within the flow field, it is necessary (as is shown in ref. 1) to allow the veloc- 
ity to be discontinuous across S. For this reason the streamline S will be called the 
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slip line. The region within the jet and orifice is denoted by D + and the remaining re- 
gion of the flow (i. e., the main stream) by D". Since the velocity (and as a conse- 
quence, the velocity potential) is discontinuous across S, it is convenient to use a super- 
script + to denote the flow quantities inside the jet (i. e. , in D + ) and a superscript - 
to denote those in the main stream. Thus, 


S(z) 


V(z) 

) 

jT(z) 


for z e D + 
for z ^ D - 


and 


fw + (z) 



for z e D + 
for z e D - 


Then and W + are holomorphic in the interior of D + and and W" are holo- 
morphic in the interior of D _ . 

A repetition of the argument given in reference 1 shows that Bernoulli's equation 
implies that 



ir(z s )i 2 



( 2 ) 


g 

at every point z of S. Since S is a common streamline to the internal and external 
flows, it is clear that 9*h'W* { z®) and ^/»W - (z^) are both constants. Moreover, the ar- 
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bitrariness in the definition of W can be partially removed by choosing these constants 
to be zero (ref. 3). Hence, 

9m W 4- (z) = 9mW~( z) = 0 for z e S (3) 

The remaining arbitrariness in W can be removed by choosing 

W + (0) = W~(0) = 0 (4) 

The conditions imposed on the velocity at infinity are (in view of the manner of nondi- 
mensionalization) 


i 

o 

for z - H 

(5) 

- 1 

for z — G 

(6) 


The remaining boundary conditions are that the normal component of the velocity vanish 
on the solid bounda.ries. These conditions are sufficient to completely determine the 
solution. They are summarized below for convenient reference. 



C(z) 



z es 


9m W + (z) =9mW (z) = 0 


9m ? + (z) = 0 ; 

z e EC 


9m z) ^ 0; 

z e HD 


9mC(z) = 0; 

z e EH 


9m i^ _ (z) = 0 ; 

z e GD 

J 


( 7 ) 


Asymptotic Expansions 

For small values of e the functions and W* can be expanded in an asymptotic 
power series in e. In view of the fact that the shape of the slip line depends on e, 

g 

these expansions imply that the coordinates of S, z , and the asymptotic jet width h 
must also be expanded in powers of e. Hence, 
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I[ 

C*= ?o + €? 1 + • 

W* = W Q + eWj + 

S S S 
z = z Q + ez 1 + . 

h = hg + eh^ + . 



( 8 ) 


s 

As pointed out in reference 1 the expansion of z does not imply that the complex 
variable z is being expanded. It is also shown in reference 1 that the coefficients in 
the expansion of are related to those in the expansion of W* by 


Cn = 


,± 

■1 


dW^ 

dz 


dW, 


dz 


( 9 ) 


The reason for omitting the superscript + or - in the zeroth-order terms of the 
first two expansions is that (as will be shown subsequently) the zeroth-order solutions 
are not discontinuous across the curve S and so there is a single function which is 
holomorphic in the entire flow field (of course, the same is true for Wq). 

The reference length l will now be chosen in such a way that 


Thus l is the zeroth-order asymptotic thickness of the jet. This is denoted symboli- 
cally by putting 



( 10 ) 
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The last expansion (eq. (8)) is then 


h — 1 + th^ + . . , 


(U) 


Zeroth -Order Solution 


When the expansions (eq. (8)) are substituted into the boundary conditions (eqs. (4) 
to (7)) and only the zeroth-order terms are retained, the following boundary conditions 
for the zeroth-order solution are obtained: First, the first boundary condition (eq. (7)) 
shows, as has already been anticipated, that the zeroth-order solution must be continu- 
ous across the slip line, and, hence, that it is characterized by functions which are 
holomorphic everywhere within the flow field. The remaining conditions show that 


W 0 (0) = 0 

A* w o( z o ) =0 


fz e EC 


9m C Q (z) = 0 

? 0 (z) - 0 
£ 0 (z) - 1 


for < 


z e HD 
z e EH 




GD 

for z — H 
for z — G 


(12) 


(13) 


The conditions (eq. (12)) merely serve to show that because of the manner in which the 
arbitrary constants have been adjusted in the complex potential the streamline emanating 
from the point D is to be taken as the zero streamline. 

Now the change in the stream function across the jet must be equal to the volume 
flow rate through the jet. Hence, if Q q denotes the dimensionless zeroth-order volume 
flow through the jet, and A^q denotes the zeroth-order change in the stream function 
across the jet, it is clear from the definition of stream function that 


Q 0 = A ^ 0 
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The last boundary condition (eq. (13)) shows that far downstream in the jet (i. e. , at the 
point C) the zeroth-order velocity goes to 1. In view of the normalization (eq. (11)) the 
asymptotic thickness of this portion of the jet must also be 1. It follows from these re- 
marks that Qq = 1. Hence, 


A^o = 1 (14) 

Now the boundary value problem posed by the boundary conditions (eq. (13)) is a 
simple free streamline problem which can be readily solved by the Helmholtz-Kirchoff 
technique. In fact, the solution to this problem has already been carried out by Ehrich 
(ref. 4). His solution, however, is somewhat inconvenient for our purposes. Thus, it 
will be necessary to use a slightly different approach in order to obtain a zeroth-order 
solution which is in a convenient form to use for calculating the higher-order terms. 

The procedure for obtaining the solution is (ref. 3), of course, to draw the region of 
flow in the hodograph plane and in the complex potential plane, and then to find the ap- 
propriate mapping of these two planes into some convenient intermediate plane (say, the 
T-plane). The shapes of these regions can readily be deduced from the boundary condi- 
tions (eq. (13)), and they are shown in figures 3 and 4 (we have put Wq = cp^ + i^/„ in 
fig. 3). The corresponding points in the various planes are designated by the same let- 
ters. The zeroth-order "slip line" is shown dashed in these figures since it does not 
correspond to a line of discontinuity and can therefore be ignored as far as obtaining the 
zeroth-order solution is concerned. The intermediate T-plane is chosen in such a way 


G 


( 0 . 0 ) 


G 


^0 


H -1 C 


E 

Figure 3. -Zeroth -order complex potential plane (W 0 -plane). 
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1 


~ v 0 



Figured - Zeroth-order hodograph (£ 0 -plane). 



Figure 5. - Intermediate plane (T-plane). 


that the region of flow maps into the upper half plane in the manner indicated in figure 5. 
We shall denote the real and imaginary parts of the variable T by £ and 77, respec- 
tively. The region of the T-plane into which the zeroth-order flow field interior to the 
jet maps is denoted by ^q, and the region of the T-plane into which the zeroth-order 
main stream maps is denoted by The dividing line between these two regions 

(which is being called for convenience the zeroth-order slip line even though no slip oc- 
durs in the zeroth-order solution) is denoted by 

Simple applications of the Schwartz- Christoffel and linear fractional transformation 
(ref. 5) show that the mappings which properly transform the W^-plane and the 
plane into the upper half T-plane in the manner-indicated in the figures are respectively 
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defined by 


and 


dWg 

dT 


1 T + 1 

7 T T 


for 77 s 0 


? 0 = T - for 77 — 0 
u T - A 


or, performing the indicated integration, 


W 0 = - (T + 1 + In T) - i for 77 > 0 

7 T 


(15). 


(16) 


Since the Wp-plane and the T-plane are essentially the same as those of refer- 
ence 1 , the results obtained therein can be used to show that the parametric equation for 
the zeroth-order slip line in the T-plane is 

1 T 2— + iTj 2_ e -i 2 for 0 < 77 < vr (17) 

tan 77 sin 77 


It follows from the first equation (9) that the points in the physical plane (fig. 2) are 
related to the points in the T-plane by 


z(T) = 



dW r 


C 0 (T) dT 


dT + Constant 


(18) 


Substituting equations (15) and (16) into this formula and using the fact, indicated in fig- 
ure 2 , that the origin of the coordinate system in the physical plane is to be at the point 
D result in 


z(T) 



(T - A)(T + 1) 

m2 


dT 


= A |t + (1 - a) I n T + — + (1 + A)1 + (A - l)i 
77 L T J 


( 19 ) 
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By definition (see figs. 2 and 5) 


z(a) = a + ib 


Hence, equation (19) shows 

a + ib = - [(1 - A)ln A + 2(1 + A)] + i(A - 1) 

7T J 


Therefore, on equating real arid imaginary parts, 

a = - [(1 - A) In A + 2(1 + A)] 


( 20 ) 


b = a - 1 


Formulation of First-Order Problem in Physical Plane 

The mapping T — Z defined by equation (19) maps the upper half T-plane approxi- 
mately into the region of flow in the physical plane. The domain <2+ is mapped into the 
crosshatched region of the physical plane shown in figure 6. The curve y Q is mapped 
into the dashed boundary S Q of this region. This region, of course, differs from the 
true interior of the jet whose boundary S is indicated by the solid line (curved) in fig- 
ure 6. 

Now the first group of the boundary conditions (eq. (7)) is specified on the curve S 
in the physical plane, whose shape is not known at this stage of the solution. As ex- 
plained in reference 1, however, these boundary conditions, correct to terms of the 
order e, can be transferred to S n by relating the values of and W* at an arbi- 

Q U g 

trary point z of S to their values at some neighboring point Zq of Sq by performing 



Figure 6. - Comparison as zeroth order and true jet boundaries in physical plane. 
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a Taylor series expansion of these quantities about 


s% s ) = s 



z 8 . Thus, 



w± (z S ) = W ± (z 8 ) + ? ±(z 8 )(z 8 - z 8 ) + . . . 

As in reference 1 substitution of the asymptotic expansion (eq. (8)) into these Taylor 

2 

series yields, after neglecting terms of 0(e ), 



W*(z S > = W 0 (z») + £ [w*( Z S) + ? 0 ( z o) z i] + 0 (£ 2 ) ( 22 ) 


And these expressions relate the values of the dependent variables W* and at the 

points of the unknown boundary S to their values on the known boundary Sq with an 
error of order e 2 . An easy calculation carried out in reference 1 shows that equa- 
tion (21) implies 


cV> 


= Ko( z o)| 2 + 2£ I<o( z o)| 2/? * 


w*?) <«(■§) 



z=z r 


+ O(e^) (23) 


Finally, substituting the expansions (eqs. (22) and (23)) into the first group of bound- 
ary conditions (eq. (7)) and equating the coefficients of e to the first power yields the 
following first-order boundary conditions on the slip line: 


R& 


? l( Z o) ? l( Z o) 1 1 

k(4) U4)\ 2 M^)i 2 

(24) 

9 m [w;(z®) + ?„($*?] = 0 

(25) 

M w l($ + «o( z oK]' 0 

(26) 
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Thus, equations (24) to (26) are the boundary conditions for the first-order solutions on 
the boundary S ’'transferred" to the. zeroth- order boundary Sq. Hence, the first-order 
boundary value problem has been transformed from one in which the shape of the bound- 
aries is unknown to one in which it is known. Notice, however, that these boundary con- 
ditions involve the variables Wj, and z^. But, is completely determined in 
terms of Wj by the second equation (9). In view of this the conditions (eqs. (24) to (26)) 
may be thought of as two boundary conditions connecting the variable with the vari- 
able across Sq (or equivalently the variable W| with the variable W~) plus an 
equation which determines once ^ are known. Thus, subtracting equation (26) 
from (25) yields 


A*wj( z S)= A*w;( z s) 


(27) 


Then in view of the second equation (9), equations (24) and (27) are the boundary condi- 
tions on S„ which connect the solution in D + with the solution ^7 in D~, and 

U c A , g 

equation (25) serves to determine z^ once ^ is known (actually z^ will be deter- 
mined in a slightly different fashion). It is shown in reference 1 that the boundary condi- 
tion (eq. (27)) can be differentiated to obtain 


9m 




= 0 


Multiplying this by i and adding it to equation (24) yield the following single (complex) 
jump condition on S°: 


g i( z o) _ h( z o) = l l (28 . 

?o( 2 o) ?o( z o) 2 l?o( z o )! 2 

The boundary conditions for the remaining (solid) boundaries are easily deduced by sub- 
stituting the first asymptotic expansion (eq. (8)) into expressions (5) and (6) and the sec- 
ond group of boundary conditions (eq. (7)) and equating the coefficients of e to the first 
power. Thus, 
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rz e EC 

9m ?j(z) = 0< 

zero 


lz e EH 

m Cj(z) = 0 

zero 

O 

t 

S 

+ y-i 

z - H 

?‘(z) - 0 

z — G 


(29) 


Solution of First-Order Boundary-Value Problem 

The boundary conditions (eqs. (28) and (29)) completely determine a boundary- value 
problem (or more precisely, two boundary- value problems connected along the curve 
S Q ) for a holomorphic function in the region of flow in the physical plane. However, 
under the change of variable z ~ T defined by equation (19) this boundary- value problem 
can be transformed into one in the upper half T-plane (fig. 5). The boundary conditions 
in the T-plane are 


?l(T) ?j(T) _ 1 x 

2 | ^ q (T) | 2 


for Tey 0 


9m ?{(£ + iO) = 0 

for £ > - 1 

9m q(£ + iO) = 0 

T-i 

1 

VI 

±SS 

u 

o 

?“(T) - 0 

for T — oo 


Cj(0) = 0 


(30) 


Clearly, the domains of definition of £ + and are and respectively. It is 

convenient to work with the sectionally analytic function © defined on the upper half T- 
plane in terms of £* by 
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II I I III I 


ii w huh 


?T(t) 

0 + (T) = -A 


for T e 


for Tec 


?n(T) 2 


?i(T) i i 

0 ( T ) = _i + A A 

M T ) 2Cq(T) 2 


It follows from the boundary conditions (eq. (30)) that 0 must satisfy 

0 + (T) - 0"(T) = T(T) for Tey°^ 

9m 0(£ + i0) = 0 for -°° < £ < +°° >• 

0(T) - 0 for T - « 


where we have put 


T(T) = A 


1 1 

UT )| 2 d(T) 


1 

1 

2? 0 (T) 

? 0 ( T ) 


C 0 (T) ? 0 (T) 


Since can be no more singular than £q, if the asymptotic expansion is to be uni- 
formly valid it follows that 0 must be bounded on the real axis. 

The function 0 can now be constructed as follows: An investigation of the behavior 
of T(T) at T = °° and T = -1 shows that it vanishes at these points like some power of 
T. Hence, the Plemelj formulas (ref. 6) show that the Cauchy integral 


f r(r) 

27ri -V r - T 
■*0 


where the integration is to be performed along in a counterclockwise direction 
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about Sq) is a sectionally analytic function which is bounded on the real axis, vanishes 
at infinity, and satisfies the jump condition (eq. (30)). However, this function is not 
necessarily real for real values of T. But this can be compensated (as shown in ref. 6) 
by adding the function 

1_ /* T(r) dr 

2iri Jcs t - T 

to equation (34). (Notice that, if f is holomorphic in the upper half plane, the function 
f defined by f(T) = f(T) is also holomorphic there and f(£) + f(£) = f(£) + f(£) is real.) 
Thus, the function 0 with the required properties is defined by 


0(T) =— - f - 1 - f for Re T > 0 (35) 

2ni Jcs t - T 2v\ Jy, r - T 
•^0 *0 

In view of the mapping T — z defined by equation (19), this completes the solution of the 
problem since equation (31) determines in terms of the known function 0. It is also 
convenient to have an expression for in terms of 0. To this end, notice that it fol- 
lows from the first and second of equations (9) that 




± _ 

1 " 


dWj 1 

dz 


dW* 

dT 


dT 

dW n 


Hence, 


dW* _ dW Q 
dT dT C 0 


Using this in equation (31) shows that 


dWj 

dT 


dW 0 . 

— 0 + (T) 

dT 



dWg 

dT 


(36) 
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dWj 

dT 


dW 0 - 1 dW 0 
0 0 (T) + A A 

dT 2 dT 


1 - 


1 

?o< T ) 


(37) 


Notice that z = 0 when T = -1 and, therefore, equation (4) implies that W^(T) both 
vanish at T = -1. Hence, integrating equations (36) and (37) between -1 and T yields 


Wj(T) - 



0 + (T) 


dWp 

dT 


dT + - W n (T) 
2 U 


W"(T) = 



^0 

dT 


dT + - 
2 


W 0 (T) 



(38) 


(39) 


Let Q denote the dimensionless volume flow through the jet. It follows from the 

o 

definition of the stream function that to within an error of 0(e ) 


Q = JW[w 0 (-l + iO) + eWj(-l + i0)] - 9m[ W q (a + i0) + eW+(A + 10)] + 0(e 2 ) 

Or using equation (14) and the fact that W|(T) and Wq(T) both vanish at T = -1 this be- 
comes 


Q = 1 - e 9m W+(a + 10) + 0(e 2 ) 

Substituting equation (38) into this expression and using equation (14) again shows 

that 


Q = 



dWg 

dT 


dT + 0(e 2 ) 


Or substituting in equation (15) this becomes 


Q = 1 + — - — 9m 

2 7 T 



0 + (T) 


AAtl dT + 0 ( e 2 ) 

T 


(40) 
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Figure 7. - Path of integration for - f A 0 + (T) dT in T-plane. 

7T T 


In view of the singularity in the denominator this integral must first be carried out 
over the path shown in figure 7 and then the limit 6-0 can be taken. Performing these 
operations yields 


- f e + (T) T + * dT = — f 0 + (£) L±. 1 d£ 

IT J. 1 T 7T £ 


i0 + (O) 


where P. V. denotes the Cauchy principal value. Because 0 + (£ + iO) is real, equa- 
tion (40) becomes 


Q = 1 + - = e0 + (O) + 0(e 2 ) (41) 

2 

Since the flow far downstream in the jet becomes uniform and since lim £“(z) = 1, it 
follows from equation (2) (see fig. 2) z— °° 

lim £ + (z) = Vl + e 
z— c 


Therefore, continuity requirements dictate that the dimensionless net volume flow 
through the jet Q be equal to 


hVl + e 

Expanding this in powers of e we find with the aid of the expansion (eq. (11)) 
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( 42 ) 


Q — (1 + chj + - • * ) ■+ — c + . . . 

= 1 + e (t + h 0 + ■ * ■ 

Hence, equating like powers of e in equations (41) and (42) shows that 

\ = © + (0) (43) 

Since equation (19) set up a one-to-one correspondence between points of the physi- 
cal plane and points of the T-plane, it is clear that equations (31), (35), (38), and (39) 
can be used to compute the first-order perturbation to the velocity and stream function 
at each point of the physical plane. 

In view of the fact that once the shape of the jet is known it is quite easy to sketch in 
the streamline patterns, the most important quantities to be obtained from the analysis 
are the shapes of the curve S in the physical plane (see fig. 2). However, since the 
viscous spreading of the jet is controlled by the pressure (or equivalently, the velocity) 
distribution along the slip line S, that quantity is also of some importance. Hence, ex- 
plicit formulas will now be obtained for these quantities by using the formulas previously 
derived. 


Computation of Boundary Values 

In view of the one-to-one nature of the mapping involved it is clear that, if 

jWw+tz) = 0, then z must be a point on the streamline which passes through the point 

D in figure 2. In addition, since the velocity potential is increasing in the direction 

D — C along S, it is clear that, if Re'W + (z) > /?gW + (0) = 0 (see eq. (4)), then z must 

be a point of the slip line S. It is clear that 9m W n (z?) = 0 and Re W n (z?) > 0 for any 
S u u u c s 

point Zq e S Q . In view of these considerations it follows that the point z* = Zq + ez^. 

will be on the first-order position of S if z^ satisfies the equation 

W*(z S ) = W 0 (z®)(l + (44) 

It is clear from this equation that when z? = 0, W + (z S ) = 0, and that Re W + (z S ) — °° as 

s s u s 

Zq — oo. Hence, the point z traverses the streamline S as Zq traverses the zeroth- 
order slip line Sq. 
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By substituting equation (44) into the expansion (eq. (22)), the first-order distance 

Q 

Zj from the zeroth- order slip line to the slip line is found to be 


— L- |wJ(T) - i W 0 (T)j for T e y Q (45) 


where the fact has been used that the curve S Q in the physical plane is the conformal 
image under the mapping T — z defined in equation (19) of the curve y Q in the T- 
plane. This also shows that 


z S = z(T) + ez® + 0(e 2 ) for T e y^ 


In addition, equation (23) shows that the magnitude of the velocity at each point z 
of the slip line is given to within terms of order e 2 by 


( 46 ) 

S 


r<z S )| 2 = i ? 0 (T) | 2 + 2e | C q (T) | 2 Re- 


?t(T) 


d? 


? 0 (T) ? q (T) dz 


°z S 

Z 1 


for T e y Q (47) 


Substituting equation (45) into (46) and substituting equation (38) into the resulting ex- 
pression yield 


„ / dW n 

z b = z(T) / 0 + (T) — - dT for T e y n (48) 

?0 (T) dT 

where the integral can be taken along the curve y Q if 0 + (T) is interpreted as the 
limiting value of 0(T) as T approaches y Q from within 

It follows from the first equation (9) and equations (46), (48), and (31) that 


gj( T > t 1 d ^o z s 

? 0 (T) + ? q (T) dz 1 


5lCT) + 5 o Z S 
C 0 (T) + dW 0 1 


= I + 0+(T) 
2 



dW r 


dT 


dT 
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Substituting this into equation (47) shows that 


|? V )| 2 = |? o < T )! 2 ( 1 + e 


2e Re 


0 + (T) - 


f\ + 

dW o J_ x 


dW 0 

(T) H dT 

dT 


for T e (49) 


The distance S measured along the curve S is given by 


S = H, 


/ 


dz 

dT 


I dT I for T e y r 


(50) 


where the integral is taken along the curve y^. Now differentiating equation (46) and 
inserting equation (45) yield 


dz 

dT 


= $ 2 - + - (wt - A w n \ 

dT r 2 dT V 1 2 u / 

^0 


C 0 \ dW 0 


dWp 

dT 


dW r 


dT 


. € Xl d ^0 1 S' 

1 + € j + Z. 

2 \? n dz Cn 


Comparing this with equation (47) shows that 



Substituting this result into equation (50) shows that 



( 51 ) 


24 



All necessary results have now been obtained. However, it is convenient to rewrite 
some of these in more explicit form. 


Explicit Formulas for Calculating Boundary Values 

Substituting equation (16) into equation (33) yields 


T(T) = -i 9m 

T T 

Or using equation (17) 

r(rj) = -i (rj + A sin e*^)A ^ — ^ ^ for 0 < 17 < n (52) 

where r(rj) is used in place of r j^-(^/sin T])e”^J . Applying the Plemelj formulas 
(ref. 6) to equation (35) shows that 

0 + (T) = - 1 - T(T) + f - _I_ f r ( T ) dr for T e 

2 27ri J</> t - T 27ri Jcs r - T 

•^0 *^0 

where the integration is to be performed in a counterclockwise direction along In 

view of equation (24), however, this can be written as 



1 

27TX 



r c iy . U c -iT7 \sin y ) 

sin y sin 77 


for 0 < T) < 77 


where, for brevity, 0 + (r;) is used in place of 


o + [(-f?/ sin v) e 


Since 
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( 53 ) 


= -ft 

\sin y / \ 


2-^os r sin y + 


sin^y 


the foregoing equation can be written as 


© + (rj)=-r( T? ) + P^- 

2 27ri 





y c 

r (y) 

1 

sin^y 


— — (y - COS y Sin y) + i 
sin 2 y 

y cot y - 77 cot 77 - i(y - 77) 


d>/ 


(y - COS y sin y) - i 


y cot y - 77 cot 77 + i(y + 77) 


dy 


for 0 < 77 < 7r ( 54 ) 


For T = 0 equation ( 35 ) becomes 


e + (0) = - A f r(y)/ i - cot y - i)dy 
Ti Sq \y } 


( 55 ) 


s s 

Upon defining z^iri) and Cq(t?) by 


-V 


z 0 (t?) = z — e 

\sin r\ 


\sm 77 / 


> for 0 < 77 < 7T 


and using equations ( 15 ), ( 17 ), and ( 53 ), equation ( 46 ) becomes 

JM. 


Z S = Zq(t]) - e 


Cnto) 


for 0 < TJ < 77 


( 56 ) 


( 57 ) 


where 
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i(t?) h r 0+ (>') J (r) d >' 


for 0 < 17 < 7 r 


(58) 


and 


J(y) s — j"(l - y cot y)^ + y^l for 0 < y < n 

Try L J 


(59) 


Upon defining VZ,(ri) by 

o 

V>) sV ooKV)| for 0 < 77 < 7 f 

and using equations (15), (16), (17), and (53), equation (49) becomes 


(60) 


v>) 

V 


= i I 2 { 1 + e + 2e [ 0+ (i?) - A(t7)I(t7)J| 


for 0 < 77 < tt (61) 


where 


a(t?) i 


7 T A sin 77 


(sin rj - 77 e lr ^)(rj e 1T ^ + A sin 77 ) 


for 0 < rj < 77 


(62) 


Let p Q be the pressure far inside the orifice (the point H in fig. 2). Since the velocity 

is zero there, it is clear that p Q = P^. Hence, it follows from equation (2) that the 

pressure at the point (x_, y_) on the slip line p(x y ) is given by 

s s s s 


p 0 - p(x s-/y 

ipv* 

2 


+ ,_S„2 = M 


= r(z u ) 


(63) 


Hence, let the pressure coefficient on the slip line Cp g be defined by 

_ P 0 -P(x s ,y s ) 

ps , o 

A pV^ 

2 P 00 


(64) 
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Upon using equations ( 15 ), ( 17 ), ( 53 ), and ( 60 ), equation ( 50 ) becomes 


S(T 7 ) _ 

H 0 V - 


(l +~\ f — ^—J(r)dy 

V z)J 0 v >) 


for 0 < 77 < 7r 


For convenience, the most important equations of this section are now summarized. 


SUMMARY OF EQUATIONS 


a = — [(1 - A) In A+ 2(1 + A)] 


> for 0 < A 


b = A - 1 




? 0 (T) = 


T - A 


for 9 m T > 0 


»[ T 


>)j + i( 


z(T) = - T + (1 - A) In T + — + (1 + A) | + i(A - 1 ) 

T 


h = 1 - — 9*n f T(ri)(- - cot 77 - Adr/ 
n \rj J 


T{ v ) - -iA(ij + A sin 77 e iT? ) 

2 

V 


0 + (T7) = i Tin) + P 


P. V. f 

2tri / 
«/ 0 


r (y)M(y)dy 


y cot y - rj COt 77 - i(y - 77) 


hri / 

‘ / 0 


_T(y)M(y)dy 


for 0 < 77 < 7r 


27 ri / y cot y - 77 cot tj + i(y + 77) 


for 0 < 77 < n 


where 



M(t 7 ) h — - — (y - cos 77 sin 7]) + i 
sin% 


( 67 ) 


zg(u) - z - -J}- e-M ^ 

\ sin 7] ) 


> for 0 < Tj < 7T 


tgOU-to- 

\ sin r\ ) 


(56) 


vj(n) = k + (z s )|v 0 


J(y) = — f(l - r cot y ) 2 + y 2 l 

7iy L J 


(60) 

(59) 


A (77) = 


2 

77 A sin 77 


(sin 77 - 77 e lr *)(rj e 1T ^ + A sin 77 ) 


for 0 < 77 < 77 (62) 


I(t?) = f V © + (y)J(y)dy 

•/0 


( 58 ) 


, s - - « -M. 

?o‘i) 


for 0 < 77 < 77 


(57) 


V^(77) 


V„ 


= C ps = ^0 


(77) | 2 jl + e + 2 e / 2 e^=> + (ri) - A (77)1(77)^ 


for 0 < r] < 7 t 


(61) 


V =o H 0 


= 1 + 


2 U v >> 


J(y)dy for 0 < 77 < 77 


(65) 


where equations (55) and (43) have been substituted into equation (11) to obtain equa- 
tion ( 66 ). 
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RESULTS AND DISCUSSION 


The numerical calculations were performed by using complex arithmetic. Hence, 
there is no need to separate the real and imaginary parts of the various formulas given 
in the preceding section. Equations (20) are used to calculate the orifice offset ratio 
B/A = b/a for various values of the parameter A. However, it is more convenient to 
present the results in terms of the orifice orientation angle defined as tan~*B/A. A 
plot of the orifice angle against the parameter a is presented in figure 8. The orifice 
angle completely fixes the geometry of the problem. Hence, once the geometry of the 
orifice is set, the parameter A can be determined from figure 8. This parameter is 
the one which appears naturally in the formulas which are used to calculate the various 
physical quantities of interest. The only other parameter appearing in the problem is e 
which gives a measure of the difference between the total pressure in the jet and the 
total pressure in the mainstream. This parameter is defined by equation (1) as 


e = 



Equation (52) is used to calculate r( 7 ?) for various values of A and these values of 
r(rj) are used together with equation (67) to calculate 0 + (£) and h for various values of 
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A from equations (54) and (66), respectively. All the physical quantities presented in 
the plots are determined by these latter two quantities. 

Now for two-dimensional jets, the jet contraction ratio is defined as the asymptotic 
jet thickness divided by the length of the orifice. Hence, the jet contraction ratio is 

JH = h 

Va 2 + B 2 Va 2 + b 2 

Substituting equations (20) and (66) into this formula gives the jet contraction ratio as a 
function of A and e or in view of figure 8 as a function of tan - *B/A and e. These 
results are presented in figure 9. It can be seen from figure 9 that, as in the case of 
the separated jet discussed in reference 1, for positive values of the orifice angle small 
changes of e result in large changes in the contraction ratio, the effect becoming more 
marked as the orifice angle is increased. The opposite conclusion holds for negative 
values of the orifice angle. However, this effect is not nearly so marked as that of ref- 
erence 1. Figure 9 also shows that, as in the case of the separated jet, for a given ori- 
fice angle increasing e always results in an increase in the jet contraction ratio. This 
increase is negligible, however, for orifice angles less than -1. 6 radians. Figure 9 
shows that the jet contraction ratio is a maximum for an orifice angle of -1.2 radians 
and falls off markedly when the orifice angle is changed. 

The parametric equations (with parameter 77) for the slip line are obtained by sub- 
stituting equations (19), (16), and (58) into equation (57) and using definitions (56) and the 



Value of 
€ 
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expression for 0^(4) discussed previously. The resulting expression determines the 
boundary of the jet. The shapes of the jet boundaries for various values of the param- 
eters e and B/A are shown in figures 10(a) to (i). Figure 10(a) corresponds to a jet 
injected normal to the mainstream B = 0. Figures 10(b) to (e) are for negative orifice 
angles (i. e. , jet injected downstream) and figures 10(f) to (i) are for positive orifice 
angles (i. e., jet injected upstream). The configurations shown in figures 10(g) to (i) 
may be strongly modified by viscous effects. Figures 10(a) and (f) show that, as in the 
case of the separated jet discussed in reference 1, when the orifice angle is greater than 
or equal to zero a small change in the total pressure within the jet results in a fairly 
large change in both the jet penetration and jet thickness. This effect becomes more 
pronounced as the orifice angle is increased. However, the effect is not as marked as 
in the case of a separated jet. The figures also show that turning the jet into the main 
stream tends to markedly decrease the flow in the jet. The extreme sensitivity of the 
flow configuration to e at large positive orifice angles indicates that the perturbation 
analysis will break down when the orifice angle is large enough. 

The pressure coefficient on the slip line is obtained as a function of the distance 
along the slip line in parametric form from equations (61) and (65) after using definitions 
(52), (54), (56), and (58). The results of these calculations are shown in figures 11(a) 
to (h). Each figure is drawn for a different orifice angle. These curves contain all the 
information necessary for calculating the viscous boundary layer along the slip line. 

The curves show that for a fixed orifice angle the velocity at both the upstream edge of 
the orifice and at the downstream end of the slip line increases with increasing e. For 
negative orifice angles the velocity tends to be relatively constant along the slip line, ex- 
hibiting a slight dip at the upstream edge of the orifice. As the orifice angle is increased 
toward zero the variation of velocity along the slip line becomes more pronounced. For 
nonnegative values of the orifice angle there is a definite peak in the velocity profiles 
which becomes more marked as the orifice angle is increased. This peak is attributed 
to the fact that the velocity is infinite at the downstream edge of the orifice. Since this 
point moves closer to the slip line as the orifice angle is increased, the velocity along 
the slip line becomes more peaked as the orifice angle increases. The contrast between 
the velocities (pressure coefficients) on the slip lines of the attached and separated jets 
is illustrated in figure 12 for zero orifice angle. The pressure coefficient of the sepa- 
rated jet rises monotonically to 1 far downstream. This smoother behavior can be at- 
tributed to the presence of the wake which adjusts in shape to keep the velocity of the 
turning jet from getting too large and the pressure from dropping too far below the static 
pressure of the stream. If the downstream wall is actually very thin, viscous effects 
will very likely cause the formation of a separation bubble at the lip. In that event, the 
velocity distribution on the slip line will fall between the two extremes shown in the fig- 
ure. 
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Dimensionless coordinate, Y/H, 


Value of 



(a) Orifice angle, tan" 1 B/A = 0. 



(b) Orifice angle, tan' 1 B/A = -2 tt/ 3. Absolute value of epsilon, I e k 0.4. 



<c) Orifice angle, tan 1 B/A = -i r/2. 



Dimensionless coordinate, X/Hq 


(d) Orifice angle, fan" 1 B/A = ~?r/3. 

Figure 10. - Jet contours for varying orifice angles and varying values of 
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CONCLU DING REMARKS 

A procedure developed in reference 1 has been applied to obtain a solution to the 
problem of a two-dimensional inviscid jet injected from an orifice at an oblique angle to 
a moving stream for the case where the jet does not separate from the downstream edge 
of the orifice. The analysis shows a peaking of the pressure coefficient along the slip 
line which increases as the orifice is tilted into the stream. Tilting the jet into the 
stream also makes the jet contraction ratio more sensitive to changes in total pressure 
in the jet. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 12, 1969, 

129-01. 
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